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ABSTRACT 

We report the results of a near-infrared survey for long-period variables in a field of view 
of 20 arcmin by 30 arcmin towards the Galactic Centre ( GC). We have dete cted 1364 vari- 
ables, of which 348 are identified with those reported in iGlass et all (1200 ll) . We present a 
catalogue and photometric measurements for the detected variables and discuss their nature. 
We also establish a method for the simultaneous estimation of distances and extinctions us- 
ing the period-luminosity relations for the JHK S bands. Our method is applicable to Miras 
with periods in the range 100 - 350 d and mean magnitudes available in two or more fil- 
ter bands. While J-band means are often unavailable for our objects because of the large 
extinction, we estimated distances and extinctions for 143 Miras whose H- and X s -band 
mean magnitudes are obtained. We find that most are located at the same distance to within 
our accuracy. Assuming that the barycentre of these Miras corresponds to the GC, we esti- 
mate its distance modulus to be 14.58 ± 0.02 (stat.) ± 0.11 (syst.) mag, corresponding to 
8.24 ± 0.08 (stat.) ± 0.42 (syst.) kpc. We have assumed the distance modulus to the LMC to 
be 18.45 mag, and the uncertainty in this quantity is included in the systematic error above. We 
also discuss the large and highly variable extinction. Its value ranges from 1 .5 mag to larger 
than 4 mag in Ak b except towards the thicker dark nebulae and it varies in a complicated 
way with th e fine of sight. We hav e identified mid-infrared counterparts in the SpitzerflRAC 
catalogue of Rami rez et alj d2008l) for most of our variables and find that they follow rather 
narrow period-luminosity relations in the 3.6 to 8.0 fim wavelength range. 

Key words: Galaxy: centre - structure - infrared: stars - ISM: extinction - stars: AGB and 
post-AGB - stars: variables: Miras 



1 INTRODUCTION 

The detailed structure of the Galaxy remains to be revealed. It is 
difficult to get a clear picture of its shape because the Sun and the 
Earth themselves are located within the Galactic disc. Furthermore, 
strong interstellar extinction towards the Galactic disc hinders 
our understanding of Galactic structure. Several kinds of tracers 
have been used to investigate this structure using vari ous observa- 
tional methods: inter stellar gas (Nakanishi and Sofue,l2003l.[2006t) 
and r e d-clump stars jNishiyama et al]|2005l : iBabusiaux & G ilmore 
120051 : Rattenbur v et alj|2007t) . for example. The tracers we make 
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use of in the present paper are Miras, i.e. long-period variables with 
large amplitudes (V > 2.5 mag) and relatively regular light curves 
(P > 100 d). 

The period-luminosity relation (PLR) of the Miras has been 
widely use d as a dist a nce in dicator since IGlass & Lloyd Evansl 
dl98lh and iFeast et al.1 Jl989h established it for the Miras in the 
Large Magellanic Cloud (LMC). Our South African - Japanese 
collaboration has recently appli ed it to a few dw arf spheroidal 
galaxies of the Loca l Group dMenzies et alj 2008, for Phoenix; 
IWhitelocke t al. 2009, for Fornax). ReikubaT HoO^ has shown that 
the PLR can be also applied to the Miras in Cen A, a peculiar ellip- 
tical galaxy beyond the Local Group. 

We can also investigat e the distribution of Miras in the G alaxy 
by making use of the PLR ( Groenewegen & Blommaert 2005; Mat- 
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sunaga, Fukushi & Nakada, 2005, here after M05). An important 
advantage of them as tracers is that we can obtain the position of 
each individual. M05 have conducted an exploratory study using 
the OGLE-II variability catalogue and the 2MASS all-sky cata- 
logue. They presented the distribution of Miras in and around the 
bulg e, restricted however to the lines of sight of the OGLE-II sur- 
vey dWozniak et al.1120021) . Nevertheless, it was shown that such a 
map can be used to study Galactic structure directly. Based on this 
method, kinematic information can easily be combined with po- 
sitional information and it is possible to locate each tracer in 6- 
dimensional phase space. This is different from the case for red- 
clump stars, for instance, because their distribution needs to be 
solved for as a group in a statistical way. Although a survey of 
Miras requires a long-term monitoring programme, we have found 
that it is a powerful tool for revealing Galactic structure. 

In this paper we present the results of our near-infrared (near- 
IR) survey of Miras towards the Galactic Centre (hereafter GC). 
iGlass et all (2001, hereafter G01) conducted a K -band (2.2 /an) 
survey over almost the same region as ours and found 409 long- 
period variables. Our survey is deeper by more than one magni- 
tude and, more importantly, includes simultaneous monitoring in 
J (1.25 fim) and H (1.63 fim) as well as K s (2.14 /im). We have 
detected 1364 long-period variables and compiled their catalogue 
including all the available photometric measurements. From the 
catalogue we select Miras and investigate Galactic structure. Our 
first step has been to establish the validity of the Mira PLR method 
for estimating both extinction and distance. In previous studies, e.g. 
M05, researchers adopted interstellar extinction values found by 
other methods in order to estimate distances of Miras based on the 
log P-K (or -mbol) relation. However, the extinction varies in a 
complicated way across the Galactic plane and is difficult to esti- 
mate in many cases. Both extinction and distance can be obtained 
simultaneously by using the PLR in two or more near-IR filters as 
is presented below. Our second step makes use of the survey data 
to explore the distribution of Miras and to estimate the distance to 
the GC, taken as the barycentre of the stellar population. We will 
also discuss interstellar extinction towards the GC region. 



2 OBSERVATION AND DATA REDUCTION 
2.1 Observations 

We used the IRSF 1.4-m telescope and the SIRIUS near-IR cam- 
era for our monitoring survey. These were established by Nagoya 
University and National Astronomical Observatory of Japan and 
are sited at the Sutherland station of the South African Astro- 
nomical Observatory. They can take images in JHK S simultane- 
ously with a field of v iew of 7.7 x 7.7 arcmin 2 (for details, see 
iNagashima et alj|l999l and lNagavama et alj|2003h . The seeing size 
is typically 1.3 arcsec and is around 1 arcsec at its best. 

We repeatedly observed 12 fields of view around the GC. 
Their central coordinates are given in Table Q] together with the 
numbers of times each one was observed, N a ^ s . A large part of the 
data was obtained in 2005 and 2006, while additional data were col- 
lected between 2001 and 2008. Each set of observations comprises 
of ten exposures for five seconds at slightly dithered positions. The 
number of the observation sets of a given field obtained per night 
ranges from one to six. We present an H-band image (or an RGB- 
composite image in the online journal) of the observed region in 
Fig.Q]based on one set of JHK B images. 



2.2 Pipeline reduction 

The raw data were first processed by using the SIRIUS pipeline 
software (Yasushi Nakajima, private communication]^. A rough 
outline of the procedures is as follows. First, a dark image was 
subtracted from every raw image in order to eliminate the effect 
of dark current. We took ten frames with the shutter of the SIRIUS 
camera closed on every observing night. These were combined into 
one dark image. Then we corrected the pixel-to-pixel variation of 
sensitivity by dividing dark-subtracted images by a flat image. The 
flat image was produced from twilight-flat exposures. Several pairs 
of images with different sky levels were collected from twilight-flat 
images taken in sequence, and the differences between two images 
of the pairs were combined to produce a flat-field image after nor- 
malization. In order to subtract the background pattern of the im- 
ages, we took exposures of a relatively sparse region with central 
coordinates 17:11:24.7, -27:27:20 (J2000.0). We took 10 dithered 
images of the blank field and combined them without adding any 
shifts to compensate for the dithering offsets. Thus we obtained 
a sky image without stellar signals but showing the background 
pattern of the images; this sky frame we subtracted from each flat- 
fielded image. The final step was to combine the 10 dithered images 
for a target field into a scientific image. Bright stars were selected in 
the sky-subtracted images, and the offsets among the images were 
calculated. The average of the images was taken after un-dithering 
them. The scientific image thus obtained is then free of spurious 
noise caused by bad pixels and cosmic ray events. 

2.3 Photometry 

Point-spread-function (PSF) fitting photometry was performed on 
the scientific images with the DAOPHOT package in IRAlQ. For 
each field one image taken under the best conditions (of weather 
and seeing) was selected from among the jV b s images for each 
filter as a reference frame. In order to standardize the magnitudes 
we compared the photometri c result of the reference frame with 
the photometric catalogue by Nishiva ma et al. (2006a). Their data 
were taken with the same instrument as ours. We used our own 
implementat ion of the Optimal Pattern Matching (OPM) algorithm 
(Tabur 2007) to match stell ar positions in ou r images and celestial 
coordinates given by|Nishiyama et al. (2006a). The magnitudes we 
obtained are on the IRSF/SIRIUS natural sy stem, whose zero-p oint 
calibration was based on a standard star in lPersson et alj dl998h . 

We also performed PSF-fitting photometry using DAOPHOT 
for the remaining (-/V b s — 1) images, and the results were com- 
pared with those from the reference frame. The comparisons were 
again made using the OPM algorithm. The positions are consistent 
among the images to within ±0.1 arcsec. We present examples of 
the magnitude comparisons in Fig. [2] The zero points of the mag- 
nitude scale for the (iVobs — 1) images were calibrated to fit with 
that of the reference based on 3-a clipped medians. The differences 
between magnitudes in the frames are distributed around zero for 
vast majority of stars, whilst the differences for variable stars fol- 
low their variations. 

1 The software is a set of scripts making use 
of IRAF (Imaging Reduction and Analysis Facility). 
http://www.z. phys.nagoya-u.ac.jp/~nakajima/sirius/software/software. html 

2 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 




Figure 1. An H-band image (or an RGB composite image in the online journal) of the observed field. North in the equatorial system 
is up, and east is to the left. Solid (or red) lines indicate I = 0° and b = 0°, while a grid with spacing of 0.1° is shown with dashed 
lines (or blue ones). The covered area is about 20 arcmin by 30 arcmin. 



2.4 Variability Detection 

Using the procedures described in the preceding section, we ob- 
tained lists of the sources detected in -/V b s images for each field 
and each filter. We expect that the faintest objects are detected only 
in some images due to variable conditions such as seeing. We in- 
cluded sources which were detected in 0.5-/V o b s or more images in 
our catalogue. Using N^ct to denote the number of detections of 
a given source, the condition can be written as TVdct ^ 0.5iV o bs- 
The number of detected sources, N Brc , is listed in TableQ]for each 
field. The total number for the 12 fields is also given. Fig.[3]shows 
histograms of JHK S magnitudes and illustrates the detection limit. 
The shaded area is the histogram for sources with TVdct 0.8-/V o b s 
while the outline is that for all the sources (iVdet ^ 0.5-/V o b s )- For 



a wide range of magnitudes, it is clear that the sources in our cat- 
alogue appear in almost all of the N h a images. The numbers of 
detections start to decrease at the faint end. The detection limit for 
each image depends on the seeing and other weather conditions as 
well as the crowdedness, which varies considerably as seen in Fig. 
[T] We here define the typical limiting magnitude as that at which 
the number of sources with iVdet 1? 0.8-/V o b s becomes less than 
half of those with iV dct > 0.5A^ obs ; J = 16.4, H = 14.5, and 
K a — 13.1 mag. For sources fainter than these magnitudes it gets 
difficult to keep the detection ratios, N^et /-^obs, l ar g er than 80 per- 
cent. The limits are indicated by the dotted lines in Fig. [3] 

We also found that the brightest sources are affected by devia- 
tions from linear response of the detectors. The deviation gets larger 
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Table 1. The observed fields. The central coordinates and the numbers of observations -/V obs are indicated. We also list the number of sources (N SIC ) 
and variables (N VS . T ) detected in each field as well as the total number in which we do not include duplicate detections in neighbouring fields. 



Field 


Centre 


(J2000.0) 


Centre (Galactic) 






^obs 










JV var 




RA 


Dec 


in 


6(°) 


2001 


2002 2004 


2005 


2006 


2007 


2008 


Total 






1745-2900A 


17:46:10.5 


-28:53:47.8 


+0.0970 


-0.0831 


1 


1 


47 


42 


2 




94 


9,266 


128 


1745-2900B 


17:45:40.0 


-28:53:47.8 


+0.0391 


+0.0119 





1 


47 


41 


2 




92 


6,797 


132 


1745-2900C 


17:45:09.5 


-28:53:47.8 


-0.0189 


+0.1068 




1 


45 


40 


2 




90 


6,041 


89 


1745-2900D 


17:46:10.5 


-29:00:28.0 


+0.0020 


-0.1409 




1 


46 


42 


2 




93 


8,692 


104 


1745-2900E 


17:45:40.0 


-29:00:28.0 


-0.0558 


-0.0460 




1 


45 


41 


2 




91 


8,451 


225 


1745-2900F 


17:45:09.5 


-29:00:28.0 


-0.1137 


+0.0488 




1 


43 


41 


2 




89 


8,185 


151 


1745-2900G 


17:46:10.5 


-29:07:07.8 


-0.0928 


-0.1987 




1 


41 


41 


2 




87 


9,576 


101 


1745-2900H 


17:45:40.0 


-29:07:07.8 


-0.1506 


-0.1039 




1 


43 


41 


2 




89 


7,399 


123 


1745-29001 


17:45:09.5 


-29:07:07.8 


-0.2085 


-0.0091 




1 


43 


36 


2 




83 


8,747 


122 


1745-2940G 


17:46:10.5 


-28:47:07.9 


+0.1919 


-0.0254 




1 1 


41 


38 


2 




85 


8,885 


174 


1745-2940H 


17:45:40.0 


-28:47:07.9 


+0.1339 


+0.0698 




1 1 


40 


39 


2 




85 


9,237 


130 


1745-29401 


17:45:09.5 


-28:47:07.9 


+0.0758 


+0.1648 




1 1 


15 


40 


2 





60 


8,636 


80 


Total 
























81,019 


1,364 
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Figure 2. Examples of magnitude comparisons between the refer- 
ence data (MJD rc f = 53665.7594 in this case) and repeated ob- 
servations. These plots are for the 1745-2900A field and the MJDs 
of the target data are indicated on top of the panels. 
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Figure 3. Histogram of JHK B magnitudes of detected sources in- 
cluding non-variables. The shaded area is the histogram for sources 
with N dct ^ 0.8N ohs while the outline indicates that for all the 
sources included in our source list (N dct ^ 0.5N o ^ s ). The solid 
lines indicate the saturation limits while the dotted line do the de- 
tection limit (see text for their definitions). 



than 5 percent for objects brighter than J = 9.0-9.5, H = 9.0- 
9.5, or K B — 8.5-9.0 mag. This limit is also dependent on seeing 
size and background brightness. In the following we do not con- 
sider magnitudes brighter than 9.5, 9.5, and 9.0 mag in J, H, and 
K B , respectively. We did not include measurements above the upper 
brightness limits in N^ct- 

Each source included in our list has magnitudes obtained for 
iVdet images in each band, where AT det can be different for each 
of the JHK S bands. We calculated the standard deviations (SDs) 
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Figure 4. The standard deviation of JV^et magnitudes is plot- 
ted against the mean magnitude for each object. The solid curves 
are the median errors at a given magnitude and the limits of non- 
variable stars are taken to be the dashed curves. The latter are 
shifted by a factor of three from the solid ones. See text for details. 



of the iVdot magnitudes, which are plotted against the mean magni- 
tudes in Fig.|4]for the case of the 1745-2900A field (as an example). 
As may be expected, the fainter sources tend to have the larger SDs. 
Assuming that the majority of the sources are not variable, at least 
to our photometric accuracy, their SDs arise from photometric er- 
rors. We estimated the error as a function of mean magnitude by 
taking the median of the SD at each magnitude. The median was 
calculated for every 0.1 mag interval, with a width of ±0.2 mag. 
This estimation was carried out from the faint to the bright end 
and we imposed the constraint that the error curve is a monotoni- 
cally increasing function of magnitude, that is, a(mi) < (j(mi) 
if mi < m2 and a (mi) = 17(7712) otherwise. Examples of the 
error curves obtained are indicated by thick solid lines in Fig. [4] 
The vertical dotted lines indicate the limit brighter than which the 
magnitudes are affected by the deviation from linear response. The 
measurements over these limits are nevertheless plotted in Fig.|4]to 
illustrate such an effect. The stars affected have large SDs, as seen 
in the bottom (K s ) panel, and are not considered further. 

The SD values above which objects were taken to be candidate 
variables were obtained by multiplying the error curve by a factor 
of three and are plotted as the thick dashed lines in Fig. [4] For ev- 
ery candidate, we inspected the light curve and the appearance of its 
image by eye in order to determine whether the variation is real or 
spurious. The candidates excluded after inspection often lie near the 
edges of the images, causing their photometry to be affected. Some 
others are severely distorted by neighbouring bright star(s) that are 
saturated. At the same time, we searched for short-period variables 
with 0.1 < P < 60 d among the candidate variables, which in- 
clude Cepheids and eclipsing binaries. We found approximately 50 
short-period variables. They are not included in this paper and will 
be reported in a forthcoming paper. The number of the remaining 
variables detected in each field, N va . T , is given in Table [T] When 
variables were found in two or more neighbouring fields near their 



edges, we confirmed that the light curves agreed and adopted only 
one entry for each variable. Thus 1364 long-period variable jf] were 
selected for our catalogue. 



3 CATALOGUE CONSTRUCTION 

In this section we discuss the estimation of periods, mean magni- 
tudes, and their accuracies. We f urther carry out cross identifica- 
tions with the surveys by G01 and lRamfrez et al. The com- 
parison with the G01 catalogue is, in some cases, useful for period 
estimation. 



3.1 Comparison with lGlass et alj J200lh 

Miras listed in G01, also see lGlass"et ai]( l2002h . were searched for 
in our catalogue. In order to adjust any systematic shifts in position 
between the two datasets, we added offsets to the positions in the 
G01 catalogue for the different G01 fields (cZra and dncc, given 
in Table |2). The application of these offsets led to better positional 
agreements. We found 348 matches with a tolerance of 3 arcsec. 
The standard deviations of the differences in Right Ascension (RA) 
and Declination (Dec) are approximately 0.5 arcsec after the offsets 
mentioned above. 

G01 listed 12 duplicated detections of the Miras in their neigh- 
bouring fields in their table 3. In course of the cross identification, 
we found several more pairs of objects listed twice by G01. Their 2- 
18 and 3-7655 corresponds to 17453194-2857473 in our catalogue, 
while 4-6/19-233 to 17453003-2905100, 9-75/12-46 to 17461431- 
2854088, and 23-1198/24-2 to 17450596-2857483. We also found 
that 2-3 is identical to 16-24/17-1. The counterpart of 16-49 is 3- 
266, not 3-226 as listed in their paper. Additionally both 9-4 and 
12-2 are close to one of the IRSF variables, 17461373-2855143, 
although they are slightly separated from each other. Their param- 
eters are close to each other and consistent with the IRSF counter- 
part. We conclude that they are also the same object detected twice 
in neighbouring fields. 

Among 56 G01 objects whose counterpart was not found, 39 
are located outside our fields of view. Coordinates of the other 
17 objects are within our survey area, but they are not detected 
as variable stars. We list those not identified in Table [3] with the 
suspected reasons for the failure. For example, the coordinate of 
6-112 is rather close to a very bright star which is saturated in our 
images, while 6-112 is listed as rather faint. Its Q flajjf] given in 
the G01 catalogue is low and their photometry must have been af- 
fected by the saturated star. Both 12-799 and 12-1236 are separated 
from 17462073-2853253 by about 3 arcsec. Different parameters 
are given in G01 for these objects, and it is not clear which one 
corresponds to the IRSF variable. Their Q flags are low. We did not 
consider the counterparts for these two G01 objects. 

We plotted a histogram of K/K s mean magnitudes which ap- 
pear in both the G01 catalogue and ours in Fig. [5] We used only 
good estimates of K B magnitudes for which Mflag = is given 
(see section [3~3t , Although we cannot directly compare the magni- 
tudes obtained with different filters, K and K a , their difference is 

3 We designate variables which are not of short period (P < 60) as long- 
period variables even if their light curves are not truly periodic. Also see 
section l4T4l regarding the definitio n of Miras. 

4 This pair was also identified by Glass et al. (2002). 

5 G01 gave a flag between (low) and 3 (high) for each object to classify 
the light curves according to their quality. 
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Table 2. The 25 fields in the G01 catalogue and the number 
of variables. Nqqi - the number of objects reported in G01; 
JV mat - the number of the IRSF ones identified with those in 
G01; dpjA an d d-Dcc - offset between the coordinates in G01 
and our catalogue. The total numbers after removing duplicate 
entries from neighbouring fields are listed in the last line. 



Table 3. G01 variables whose counterparts were not found in 
our survey. Superscripts (a-d) indicate our suggested explana- 
tions. 



Field Unmatched objects 



3 


2389 d , 2832°, 2834° 




5 


59°, 164 a ,2856 a 




6 


112 d 




7 


361 d 




8 


5 a ,31 a ,97 b 




10 


60 b 




11 


23 < \27 a ,2449 a ,4503 i > 




12 


2 d , 6 a , ll a , 136 a , 228 a 


, 799 d , 1236 d 


13 


30 a ,73 a 




14 


150 a ,463 a 




15 


5 a , 10 a , 26 a 




16 


2993 d 




17 


118 c 




19 


685 d 




20 


ll a , 22 a , 46 a , 70 a , 99" 


, 116 a , 133 6 , 


21 


6°, 8°, 17 a , 38 a , 185 a , 


5719 a 


22 


31 d ,60 a , 100 a 




23 


5 a , 15 c ,42 a ,50 a , 114 a 




24 


29 a 




25 


7" 





Field 


JVgoi 


JVmat 


cZra (arcsec) 


(foec (arcsec) 


GC1 


14 


14 


—0.57 


1.01 


GC2 


26 


26 


—0.35 


1.16 


GC3 


35 


32 


0.39 


0.50 


GC4 


13 


13 


0.45 


1.55 


GC5 


10 


7 


—0.30 


1.14 


GC6 


21 


20 


— 0.65 


0.89 


GC7 


15 


14 


-0.54 


0.74 


GC8 


12 


9 


-0.39 


0.60 


GC9 


23 




— 0.99 


0.97 


GC10 


19 


18 


-0.67 


0.51 


GC11 


11 


7 


2.64 


0.83 


GC12 


22 


15 


-0.31 


0.07 


GC13 


20 


18 


-0.34 


0.71 


GC14 


18 


16 


-0.29 


0.93 


GC15 


8 


5 


0.13 


0.78 


GC16 


26 


25 


-0.52 


0.69 


GC17 


17 


16 


-0.16 


0.81 


GC18 


7 


7 


0.13 


0.79 


GC19 


21 


20 


-0.09 


0.93 


GC20 


19 


11 


-0.17 


0.70 


GC21 


10 


4 


-0.16 


0.88 


GC22 


23 


20 


-0.20 


0.92 


GC23 


21 


16 


-0.34 


0.57 


GC24 


5 


4 


0.13 


1.25 


GC25 


6 


5 


-0.26 


0.94 


Total* 


405 


348 







* Duplicate detections in neighbouring fields are not counted. 



not expected to be large. It is obvious that many bright sources are 
saturated in our data and that many fainter ones were not detected 
in G01. With our better image quality, our detection limit is deeper 
and we also discovered many new variables which are brighter than 
the faintest G01 objects (K ~ 12), as well as fainter variables. 
Whilst many G01 variables were saturated in our A' s -band images, 
a large fraction of them were detected in the J- and/or //-bands. 
There are 348 matches between the two catalogues and none of 
these were saturated in our J-band images. It should also be men- 
tioned that saturation was not the reason for the 56 unsuccessful 
matches (see TableO. 

3.2 Estimation of periods 

We used both our own data and the G01 catalogue in order to es- 
timate the periods of the detected variables. First we carried out a 
period search of our monitoring data using the least-squares method 
to fit sinusoidal curves. We searched for periods between 100 and 
800 d fitting sinusoidal curves with minimum residuals. Then we 
examined the light curve folded with the obtained period by eye. 
We confirmed the periods for 483 variable stars. The periods ob- 
tained by us are written as Pirsf- 

There are 246 sources for which G01 obtained periods, Pgoi, 
among the variables that have Pirsf values. We have plotted 
log(P G oi/PiRSF) values against log(PiRSF) in Fig. [6] The peri- 



a The object is located out of our survey area. b A probable 
counterpart exists in our image, but we did not detect vari- 
ability. c Two or more stars are visible with approximately the 
expected brightness, but we did not detect any corresponding 
variation. d We found no clear counterpart at around the given 
coordinate. 




Magnitude {K/K s ) 

Figure 5. Histogram of K B magnitudes which appear in our cat- 
alogue (the thick outline) and K magnitudes of the G01 variables 
(hatched region) 



ods agree with each other in most cases. The standard deviation of 
log(PGoi/PiRSF) around zero is 0.025, where we have included 
those within ±0.1 dex. This indicates that the accuracies of the pe- 
riods estimated by us and by G01 are within 0.025 dex in most 
cases. We conclude that Pgoi and Prsf are consistent with each 
other when their difference is less than 0.075 dex; this level is in- 
dicated by the horizontal dashed lines in Fig. [6] When the values 
are consistent we take the period as Prsf in our catalogue. For the 
25 objects where Pgoi is not consistent with Pirsf we also use 
Pirsf- For the objects with counterparts in G01, we can combine 
their photometry with ours to estimate the period. We fitted sinu- 
soidal light curves for both A'-band (G01) and if s -band (ours) pho- 
tometry; means and amplitudes are independently obtained for two 
datasets allowing differences in K- and .Kg -band, but common pe- 
riods and dates of maxima are used for the fitting. These procedures 
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Figure 6. Comparison of Pirsf and Pgoi- The 
log (Pgoi/PrSf) values are plotted against log (Pirsf) 
whenever possible. The dashed lines indicate the range of 
±0.075 dex within which we conclude Pirsf an d P301 agree 
with each other. 



Table 4. Definition of Pflag. Number of variables with each 
Pflag is also given. 



Pflag JVvar Note 






221 


Pgoi an d Pirsf agree with each other to within 
0.075 dex. We adopt Pir Sf . 


1 


25 


Pgoi an d Pirsf disagree and we adopt Pirsf- 


2 


237 


Pgoi is not given and we adopt Pirsf- 


3 


66 


Pgoi nt our light curve but Pirsf is not prop- 
erly obtained. We adopt Pgoi- 


4 


30 


Pqoi d° es not fit the IRSF light curve. No peri- 
odicity is found. 


5 


785 


Pgoi i s n °t given and no clear periodicity is 
found with the IRSF data. 



allow us to carry out better period estimations using the photometry 
with a long time coverage of 1994-2008. We confirmed that thus 
obtained periods agree well with Pirsf to within ±0.01 dex, and 
we replaced the Pirsf whenever possible. 

There are 96 variables with Pgoi among the variables whose 
periods were not obtained by us. For 66 of them, we have confirmed 
that Pgoi can reasonably be used to explain the variations we de- 
tected. We have adopted those Pgoi in our catalogue. However, 
Pgoi does not fit the light curves for the remaining 30 variables. We 
use Pflag as defined in Table|4]to indicate the status of the periods. 
The numbers of variables with each Pflag are also given in Table 
[4] In summary, we obtained periods for 549 variables (Pflag = 
to 3). Light curves folded with the obtained periods are included in 
the online material. Fig. [7j shows examples of the light curves for 
the twenty-five cases having the shortest periods. 

3.3 Estimation of mean magnitudes 

In this paper, for simplicity, we define our mean magnitude as the 
intensity mean of the maximum and minimum of our measurements 
of a given star. We can apply this definition, unlike Fourier means, 
regardless of whether or not the light curve is periodic. Another 
possible definition is an intensity mean of all the measurements, 
but this can be affected by uneven sampling. 

In order to trust any photometric distance scale, it is important 
to obtain magnitudes for both targets and calibrators in the same 
manner. In the case of the PLR for pulsating stars, the definition 
of mean magnitude can have a significant effect. For example, in- 



tensity means differ from magnitude means. In our case, we will 
estimate mean magnitudes in the same manner as llta et al.l J2004bl. 
hereafter 104) because we will use their catalogue to calibrate our 
PLR (see section [5~TT >. 104 combined images at ten arbitrary epochs 
before carrying out photometry, so that their magnitudes corre- 
spond to intensity means of ten points. Because we have obser- 
vations at more epochs, approximately 90 in most cases, we should 
be able to get better estimates of mean magnitudes. Nevertheless, it 
is necessary to be careful about any systematic difference between 
our means and those in 104. 

We have performed Monte-Carlo simulations to assess the un- 
certainties of the ten-epoch means and their systematic difference 
from ours, if any. First, we randomly produced ten values between 
and 1 as test phases, and for each of them we obtained a mag- 
nitude from an interpolation between the closest points on both 
sides of our folded light curve. Then, we took an intensity mean 
of the magnitudes at the selected ten phases, which gives us a sim- 
ulated value for a ten-epoch mean. We repeated these procedures 
ten thousand times for each light curve to get simulated magni- 
tudes, m s i m . Then we examined how the ten-epoch intensity mean 
varies, depending on the randomly selected phases. We illustrate 
three examples of the simulations in Fig. [8] The histograms indi- 
cate how the simulated means, H s i m , are distributed around the 
means of maximum and minimum, which are indicated by the ver- 
tical dashed lines. The dispersions er s i m of the histograms suggest 
that a ten-epoch mean may have a significant uncertainty as large as 
±0.1 mag (e.g. the middle panel of Fig.[8}. Furthermore, the means 
of the simulated magnitudes (m s i m ) do not always agree with the 
vertical lines. Such differences depend on the shapes of the light 
curves. 

We now examine whether the difference between 104 and our 
results introduces a systematic offset or not. We consider only light 
curves which have good coverage in phase for the Monte-Carlo 
simulations. We also exclude light curves whose periods are not 
confirmed with our data (i.e. Pflag in section [3T2] should be 0, 1, 
or 2). There are 48 (J), 84 (H), and 85 (K B ) such light curves 
without any gap larger than 0.2 in phase among those with P = 
100-350 d. We calculate differences between the means from our 
observation and those from the simulations, i.e. m — (m s i m ), which 
are plotted against amplitudes in Fig. [9] The uncertainties cr B lm of 
a ten-epoch mean inferred from the simulations are indicated by 
error bars. One finds that both \rn — (m s i m )| and <j s i m tend to be 
large for light curves with large amplitude. The amplitudes of the 
examples presented in Fig.[8]are moderate among our sample. The 
weighted mean of the differences is close to zero. This suggests that 
the difference between our means and those in 104 is small, within 
±0.01 mag or less. 

We calculate the mean magnitudes, defined as intensity means 
of maximum and minimum, in JHK B whenever possible. How- 
ever, we did not always detect variables in all the three bands. Some 
variables are detected in H and A' s but are too faint in J, while 
some others are detected in J and H but are too bright and satu- 
rated in A" s -band. We define the Mflag as described in Table[5]in 
order to give an idea of the quality of the magnitude or a reason 
why there was no detection in a particular band. In some cases, the 
number of detections is enough to include the object in our list, i.e. 
Ndct ^ 0.5iV o b s , but a few measurements are missing at around 
maximum or minimum for instrumental reasons. We checked light 
curves for which we failed to get one or more points. The phases of 
the failures were indicated in the plots and we examined whether 
the unavailable points may have had a significant effect on the mean 
magnitudes, based on the shapes of the light curves in a given filter 
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Figure 7. Examples of periodic light curves are plotted for 25 Miras with the shortest periods. In each panel we list ID, period (d), 
and name of filter used for the plot. Plots for all of those with given period are available as online material. 
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Figure 8. Thi'ee examples of Monte-Carlo simulations of mean magnitudes from ten-epoch data are shown. For each .ff-band light curve displayed in the 
small frame, intensity means were repeatedly calculated for randomly defined groups of ten epochs. The means obtained are plotted as a histogram, while the 
thick dashed line indicates the intensity mean of maximum and minimum, as listed in the catalogue. 
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Figure 9. Expected differences between means obtained from 
observation, m, and those from the Monte-Carlo simulations, 
(m s j m ). The size of an error bar indicates the uncertainty in 

<m sim ). 



and in other filter(s) when available. Then we determined whether 
the light curves are affected (Mflag = 4-6) or not (Mflag = 0). 
The Mflag for each variable has three digits between and 7, each 
of which denotes the state of (un)available magnitudes in the order 
of J, H, and K s . 



3.4 Comparison with lRamfrez et alj J200S|> 

Ramir ez et al. N2008h compiled a catalogue of 1,065,565 sources in 
the 2.0° x 1.4° region around the GC based on their SpitzerflRAC 
survey. About 60,000 objects lie within the field we observed. They 
obtained magnitudes in four wave-passbands, i.e. [3.6] , [4.5], [5.8], 
and [8.0], so that we can investigate features of our objects at longer 
wavelengths by making use of their results. We identified 1,233 
counterparts in their catalogue with a tolerance of 0.5 arcsec. The 
coordinates in the two catalogues agree very well, with standard 
deviations of 0. 13 arcsec in RA and Dec. 



3.5 The catalogue 

We present the first ten lines of our catalogue in Table [6] as an 
example. The entire list is available in machine-readable form in 
the online version of this article. We include ID, mean magnitudes 
(JHKs), peak-to-peak amplitudes (A J, AH, AK a ), and Mflag 
(see Table |5). The magnitudes and the amplitudes are listed as 
99.99 and 0.00 respectively when we did not detect the object 
in a given band. We also list the IDs of the sources in G01 and 
lRanurezetai] lf2008) whenever counterparts were found. For those 
with period available, P and Pflag (see Table|4]l are also indicated. 
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Table 5. Definition of Mflag. The number of variables with each Mfiag is also given for J, H, and K B magnitudes respectively. 



Mflag 


Nj 


N H 




Note 





361 


990 


936 


The magnitude was obtained without any clear problems. 


1 


9 





2 


We did not obtain the magnitude because the object often falls out of the field of view. 


2 





18 


222 


We did not obtain the magnitude because the object is too bright. 


3 


833 


178 


4 


We did not obtain the magnitude because the object is too faint. 


4 


15 


35 


21 


We obtained the mean magnitude, but it is affected by its location around the edge of the images. 


5 


2 


19 


116 


We obtained the mean magnitude, but it is affected by saturation or deviation from linear response. 


6 


144 


124 


63 


We obtained the mean magnitude, but it is affected by the fact that the object was not detected in 










some cases because of its faintness. 



Table 7. The first 10 lines of the released table of light variation. This 
is a sample of the full version (1 16887 lines), which will be available in 
the online version of this journal. 



No. 


MJD 


J 


H 


K B 


0001 


52079.9757 


99.99 


11.62 


9.52 


0001 


52343.1326 


99.99 


11.62 


9.68 


0001 


53482.0590 


99.99 


11.76 


9.54 


0001 


53537.1116 


99.99 


11.79 


9.55 


0001 


53540.8221 


99.99 


11.83 


9.58 


0001 


53545.8849 


99.99 


11.84 


9.64 


0001 


53545.9694 


99.99 


11.85 


9.57 


0001 


53548.8259 


99.99 


11.81 


9.54 


0001 


53548.9591 


99.99 


11.84 


9.56 


0001 


53549.8893 


99.99 


11.84 


9.62 



Table 8. Combinations of bands in which we detected vari- 
ables. The number and percentage are given for each combina- 
tion. The values in parentheses count only variables with pho- 
tometry flags of for the given filters. 



Band combination 


Number 


Percenta 


ge (%) 


J + H + K B 


303 


(119) 


22.2 


(8.7) 


+H + K S 


655 


(485) 


48.0 


(35.6) 


J+ +K B 





(0) 


0.0 


(0.0) 


J + H+ 


201 


(152) 


14.7 


(11.1) 


J 


18 


(14) 


1.3 


(1.0) 


H 


9 


(7) 


0.7 


(0.5) 


K B 


178 


(121) 


13.0 


(8.9) 


Total 


1364 




100.0 





We also present a list of photometric measurements. The JHK B 
magnitudes obtained for all the catalogued variables are listed in 
one text file and each line includes the ID number, the modified Ju- 
lian date (MJD), and the JHK S for a given star. Table|7]shows the 
first ten lines as a sample. 



4 PROPERTIES OF CATALOGUED VARIABLES 
4.1 Magnitude and photometric flag 

Among the 1364 catalogued variables, we detected 522 variables 
in J, with 1168 in H and 1136 in K B . If we restrict ourselves to 
those with Mflag = 0, the numbers of variables are 361, 990, and 
936 in J, H, and K B respectively (Table[5j. Table [8] indicates the 
number and percentage for each combination of bands in which 
variables are detected. For example, we detected 303 variables in 
all of the JHK B filters while 178 were detected only in K B . If 
we consider only magnitudes with Mflag = 0, the above numbers 
reduce to 119 and 121, respectively. Approximately 70 percent of 
the catalogued variables have both H and A" s magnitudes (22.2 
percent with JHK B and 48.0 percent with HK B ). 

We have plotted histograms of J, H, K B , and H — K B in 
Fig. [10] Only magnitudes with Mflag = are considered. The 
hatched regions correspond to the histograms for the variables 
whose JHK B magnitudes are all available. Those with JHK B 
magnitudes occupy different ranges in the three panels. For J- and 
i/-band magnitudes, the brightest objects tend to lack A' s -band 
magnitudes (see the regions of J < 14 and H < 11). This is 
because the corresponding objects are too bright in K B . In contrast, 
J-band magnitudes are not available for a large number of objects 
which are faint in H and K B (H > 12.5 and K B > 10.5). This 
is because these objects are too faint in J, predominantly because 



they suffer from strong interstellar extinction. The bottom panel 
shows that only relatively blue variables have all the JHK B mag- 
nitudes, supporting the conclusions above. 

4.2 Colour-magnitude and colour-colour diagram 

In Fig. [TJJ we present colour-magnitude diagrams for our cata- 
logued variables (black points for those with Mflag = and black 
crosses for those with Mflag = 4-6) and non- variables (grey dots). 
The arrows indicate the reddening vector directions, with lengths 
corresponding to Ak s = 1 mag. It is apparent that the objects with 
strong extinction are often fainter than the detection limit, espe- 
cially in the J-band. 

We also plot the colour-colour diagram in Fig. [12] Due to 
large and variable extinction, the distribution of sources is elon- 
gated along the direction of the reddening vector. There is a group 
of stars which scatter towards the right side of the main diagonal 
sequence. We inspected these stars on the images and found that 
tiny displacements of centroid are recognizable among the JHK B 
images for at least some of them. In Fig. [13] we plot histograms 
of the differences between positions, both RA and Dec, obtained 
in J and K B images for individual sources. The panels (a) and (b) 
include the objects belonging to the main group and the distribu- 
tions resemble Gaussian functions centred on zero as expected. On 
the other hand, the panels (c) and (d) include objects with pecu- 
liar colours, i.e. J — H < 1.767(H — K B ) — 1.0, and show to- 
tally different distributions from those in (a) and (b). There is no 
sharp peak at zero and the differences distribute rather uniformly. 
This clearly suggests that different stars having almost the same 
positions were incorrectly combined to produce sources with pe- 
culiar colours. In other words, different stars may contribute to the 
fluxes we measured in the images in different filters: a bluer star to 
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Table 6. Catalogue of variable stars. After the numbering ID between 1 and 1364, the ID from RA and Dec (J2000.0) follows. JHK S mean 
magnitudes and peak-to-peak amplitudes are listed together with Mflag (section [33} . When the magnitudes and the amplitudes are unavailable, 
we put 99.99 and 0.00, respectively. Counterparts of Ramirez et al. 12008) and Glass et al. 12001) are indicated whenever identified. Also listed are 
periods (000 when not available) and Pflag(section l3.2t . This is the first ten lines of the full catalogue which will be available in the online version. 



No. ID Mean magnitude Amplitude Mflag Counterpart Period Pflag 

J H K B A J AH AK S R08 G01 



0001 


17445379-2857241 


99.99 


11.71 


9.51 


0.00 


0.34 


0.38 


100 


0401463 




000 


5 


0002 


17445418-2906303 


99.99 


12.50 


9.71 


0.00 


0.94 


0.71 


300 


0402528 


22-166 


339 





0003 


17445435-2857415 


99.99 


11.39 


9.15 


0.00 


0.88 


0.33 


105 


0402897 




424 


2 


0004 


17445435-2904563 


12.76 


10.05 


99.99 


0.87 


1.12 


0.00 


002 


0402921 


22-14 


280 


3 


0005 


17445442-2856527 


99.99 


12.91 


10.64 


0.00 


0.30 


0.22 


300 


0403085 




000 


5 


0006 


17445452-2856168 


99.99 


13.35 


11.09 


0.00 


0.62 


0.44 


300 


0403313 




000 


5 


0007 


17445460-2852042 


99.99 


13.87 


10.64 


0.00 


0.64 


1.12 


360 


0403552 




000 


5 


0008 


17445463-2903587 


99.99 


13.50 


10.13 


0.00 


1.68 


1.48 


360 


0403617 




457 


2 


0009 


17445482-2903208 


14.12 


10.82 


99.99 


0.53 


0.36 


0.00 


002 


0404128 




000 


5 


0010 


17445482-2905486 


16.13 


11.87 


9.46 


0.58 


0.54 


0.60 


000 


04041 17 
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Figure 10. Histograms of obtained magnitudes, J, H and K B from 
top to bottom and also that of the H — K B colour. Those for the 
objects with three-band magnitudes are shaded. 



the J-band flux and a redder star to 7\s-band. If we combine the 
flux of a foreground star on the main sequence and that of a red- 
dened giant, for example, the combined source will be located on 
the right side of the main distribution in Fig. Q2] The objects with 
J — H < 1.767(H — K s ) — 1.0 constitute not more than 3 percent 
of the number of sources with all three JHK S magnitudes, so that 
they hardly affect any of our results. 



4.3 Period distribution 

We have plotted period histograms for the variables in our cata- 
logue and also in a few others (Fig. 114b. Compared with the sample 
in G01 (the shaded area) in the top panel, we have added signifi- 
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Figure 11. Colour-magnitude diagrams, (J — H)-H in the top 
panel and (H — K B )-K B in the bottom. The black points show the 
distribution of the catalogued variables with Mflag = while the 
black crosses include those with Mflag = 4-6. Grey dots show all 
the detected non-variable sources. The arrow indicates the redden- 
ing vector corresponding to Ax e = 1 mag. 
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Figure 12. Colour-colour diagram, (H — K B ) - (J — H). Symbols 
are the same as in Fig.QT] The arrow indicates the reddening vector 
corresponding to j4#. = 1 mag. 
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Figure 13. Histograms of the difference between positions in J 
and K B images. The upper panels (a) and (b) are for objects with 
normal locations in the colour-colour diagram, while the lower pan- 
els (c) and (d) are for objects with peculiar locations, i.e. J — H < 
1.767(ff - K a ) - 1. 



cantly more variables with periods (see the thick line). The other 
panels show the histograms for the variables in Whitelock, Feast 
& Catchpole (1991, middle) and in M05 (botto m). Their objects 
are lo cated away from the GC. The sample in IWhitelock et al.1 
dl99ll) was selected from the IRAS point-source catalogue, while 
M05 i nvestigated variables found in the OGLE variability cata- 
logue (Wozniak et al. 2002). Because of the selection bias caused 
by the different wavelengths used in the surveys, the period dis- 
tributions in these two early studies are different from each other 
as discussed in M05. The former has variables with P > 600 d, 
while the latter has a more significant tail towards the short period 
end. In contrast, the histogram for our catalogue includes both the 
long-period variables and the short-period tail. It is suggested that 
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Figure 14. Histograms of periods. In the top panel, the thick line 
shows the histogram of the periods adopted in our catalogue, while 
the shaded area indicates that of Pqoi for the G01 variables with 
IRSF counterpart. In the middle and bottom pan e ls sho wn are the 
period distributions reported in IWhitelock etai] jl99l|) and M05 
respectively. 



our catalogue contains a more comprehensive sample of Miras than 
previous ones. 



4.4 Amplitude distribution 

Fig. 1151 shows histograms of peak-to-peak amplitudes in JHK S . 
The hatched region corresponds to the variables with periods (191 
in J, 432 in H and 319 in K s ), while the outline includes all the 
variables with Mflag = detected in each filter. The histogram 
of A J shows separated peaks: one for small-amplitude variables 
(A J < 0.5) and the other(s) for large-amplitude variables. Such a 
distribution is also seen in the histogram of o ptical amplitudes (e.g. 
see fig. 4 and 6 in iGlass & Llovd Evansl l2003). Many of the large- 
amplitude variables have been found to be periodic as expected for 
Miras, whilst we failed to obtain periods for many of the small- 
amplitude ones. The separation gets less conspicuous as the wave- 
length increases. The hatched region i n the AK S panel r esembles 
the distribution for Miras reported by iGlass et al, I J 19951 : see their 
fig. 3). 

It is traditional to separate the Miras from the other long- 
period variables such as semi-regulars by their large visible ampli- 
tudes, i.e. AV ^ 2.5 mag. As we will see in section |5TI we select 
periodic variables with Al > 1 mag among the LMC variables 
reported in 104 to calibrate the PLR. However, we do not have V- 
or /-band amplitudes for our objects. Most Miras have amplitudes 
larger than 0.4 mag in K (e.g. see Whitelock, Marang & Feast, 
2000), but the threshold in the near-IR amplitudes is not well estab- 
lished. It is even unclear if one can pick up the same group of Mi- 
ras by using observables at different wavelengths. There are some 
periodic variables with relatively small amplitude. For example, 
17451726-2844420 (P = 133) and 17462476-2902241 (P = 177) 
have rather small amplitudes (AJ < 0.3 mag), but their light 
curves show reasonably periodic variations (Fig. [T|>. These objects 
may be SRa type, which can be as periodic as Miras but are dis- 
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Figure 15. Histogram of amplitudes: AJ, AH and AK B from 
the top panel to the bottom. The outline includes all the variables, 
while the shaded area indicates those with period. 



tinguished by having relatively small amplitudes. However, some 
recent studies have raised questions on the distinction between the 
Miras and at least some of the semi-regulars. Lebzelter, Schultheis 
& Melchior (2002), for example, showed that a sizeable number of 
small-amplitude variables are as regular as Miras. 

The most important requirement in our selection is that our 
distance indicators follow the PLR of Miras. There are several se- 
quences in the well-known period-luminosity di agram of pulsat- 
ing red giants JWood et al.|[l998l ; llta et al.f 2004a ) . Nevertheless, it 
is known that Miras if properly chosen are exclusively found on 
the sequence C, which makes them valuable distance indicators. 
Unfortunately, it is currently impossible to tell whi ch sequence 
a given small-amplitude variable is located on as iFeasd ( 2004) 
pointed out. Some small-amplitude variables coexist with Miras on 
the sequence C. We here consider a periodic variable to be a Mira 
if its AJ, AH, or AK S , whichever may be available, is not less 
than 0.4 mag. Among our objects with periods, the number of such 
small-amplitude variables excluded is not large, about 10 percent. 
As will be discussed in section 15.2.1 1 at least some of the small- 
amplitude variables seem to be on a different sequence from that of 
Miras in the P-L diagram. 

Fig.[l6]shows i\" s -band amplitude plotted against period. One 
finds that the amplitudes do not exceed 1.5 mag for those with P = 
100-350 d, which we will use as distance indicators. This property 
does not change in the cases of A J and AH. 




2.0 2.2 2.4 2.1 

log (P) 

Figure 16. Amplitudes in the Kg-band plotted against periods. 
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Figure 17. The K B — [8.0] colours are plotted against the periods 
(top) and the amplitudes (bottom) for periodic variables. 



4.5 Mid-IR properties 

Using the counterparts identified in the catalogue bv lRarmrez et al.l 
(2008), we can study the mid-IR properties of our variables. We 
have plotted K s — [8.0] colour against period and amplitude for 
282 Miras with Mflag = in K B . In the top panel of Fig. [IT] the 
mid-IR colour becomes large at the long-period end, P > 350 d, 
similar to the amplitude trend illustrated in Fig. [16] The lower 
panel of Fig. [T7] shows that the variables with large enhancement 
in K s — [8.0] are, in fact, associated with large amplitudes. We also 
examined similar plots using [3.6], [4.5] and [5.8] magnitudes in- 
stead of [8.0], and found similar trends. A large excess in the mid 
IR is expected for a thick circumstellar dust shell produced through 
mas s-loss. This behaviour is well known from early observations 
(e.g. IWhitelock et ai] |l994l) . Detailed study of the mid-IR proper- 
ties and the mass-loss process is outside the scope of this paper. 
Please note that no correction has been made for interstellar ex- 
tinctio n and also that the mid-IR magnitudes from iRarmrez et all 
( 2008) were obtained from a single-epoch survey made with Spitzer 
Space telescope. 



Miras towards the Galactic centre 



13 



5 DISTANCE AND REDDENING 



5.1 Relation from the LMC 



In this paper we make use of the PLR for Miras in the LMC in order 
to make distance estimations. 104 published a near-IR catalogue 
of LMC variables which were detected in the OGLE-II survey by 
Udalski, Kubiak & Szymanski (1997). We select 134 Miras under 
the following criteria, 

• AI > 1, 

• J-K s < 1.4, 

• K s < 13, 

• P < 350, 

as calibrators. The colour criterion is used to exclude carbon-rich 
stars. Miras in the Galactic bulge are predominantly oxygen-rich 
(see Schultheis, Glass & Cioni, 2004, for example). The last crite- 
rion is used to exclude objects with thick circumstellar dust and hot- 
botto m-burning stars (Whi telock et alj|2003l : lGlass & Lloyd Evansl 
l2003h . The near-IR observations for the 104 catalogue also made 
use of IRSF/SIRIUS. Although 104 converted their magnitudes into 
the system of the Las Campanas Observatory (LCO), we have con- 
verted them back into the natural system of the IRSF/SIRIUS based 
on the transformation they used. Fig. [Tg] plots the JHK S magni- 
tudes against periods. The period distribution of the LMC calibra- 
tors is rather uniform except the longest period range (log P > 
2.44). In this range more c arbon-rich variables are found in the 
LMC than oxygen-rich ones (Groenewegen 2004). This situation is 
different from that of the Galactic bulge. Nevertheless, the effect 
of the period distribution in determining the PLR is small if any. 
There has been no report to offer evidence of a change in the slope 
of the PLR. 

The PLR can be determined to be: 



Mj = -2.875(logP - 2.3) - 5.798, 
M H = -3.186(logP - 2.3) - 6.508, 
M Kb = -3.555(logP - 2.3) - 6.883, 



(1) 

(2) 
(3) 



with residual standard deviations of 0.19, 0.19 and 0.17 mag, re- 
spectively. The distance of the LMC is assumed to be 18.45 mag 
(see section 1531 for its uncertainty) and extinction corrections of 
Aj = 0.06, A H = 0.03 and A Ks = 0.02 have been applied. The 
residual scatter is least in K s , but the differences are small. This 
indicates that the PLR in J and H can also be useful for distance 
estimation. More importantly, one can estimate the amount of inter- 
stellar extinction together with the distance by using the PLR in two 
or more wavebands. The residual scatter is slightly larger than that 
in iFeast et alj ( 1989) for the log P-K s relation while those around 
the PLRs in J and H are similar. The relatively large scatter is pre- 
sumably because the magnitudes in 104 were obtained as intensity 
means based on ten-epoch monitorings. The simulation presented 
in section |331 suggests that the ten-epoch mean can have an uncer- 
tainty of ±0.1 mag compared to a well determin ed value. One can 
expect that the scatter reduces to that reported in lFeast et alj (1989) 
after mean magnitudes are properly obtained. 

It is necessary to call attention to possible differences be- 
tween Miras in the LMC and those around the GC. It is suspected 
that the PLR of Miras may be dependen t on stellar parameters 
such as metallicity. Readers are referred to lWhitelock et alj ( 2008) 
for a review of the PLRs in different environments. On the other 
hand, a difference in photometric systems can also have a large ef- 
fect. Fig. [19] plo ts period-c olour relation (PCR) for the Miras in 
the LMC dFeast et alj[l989l : 104) and those in the Galactic bulge 




Figure 18. Period-luminosity relation for Miras in the LMC, plot- 
ted for 134 O-rich Miras selected from the catalogue of 104 (see 
text for the selection criteria). 



dGlass et alj [ 19951 : M05). For the relation of lGlass et aj] dl995h . 
we calculated the PCR by using the samples with period shorter 
than 350 d after two objects with extreme colours were omitted. 
Lines are drawn for a period rang e of samples in each paper; the 
shortest period in lGlass et al] d 19951) is around log P = 2.2, while 
there are several o bjec ts with period clos e to 100 d in 104 and M05. 
IFeast et alj dl989h and lGlass etall d 1995b used the Mk III photome- 
ter attached to the 1.9-m telescope at SAAO Sutherland to observe 
the Miras in the LMC and those in the Galactic bulge towards 
the Sgr I window respectively. Observations by 104 were con- 
ducted by using the IRSF/SIRIUS also at Sutherland, while M05 
used J HK S magnitudes take n from the 2MASS point-source cat- 
alogue (Skrutskie et al. 2006). Although transformations between 
these systems are available for normal types of stars, it is highly 
uncertain whether they are applicable to Miras which have strong 
absorption due to H2O and other molecules. Based on Fig. [19] it is 
difficult to discuss any population effect that may exist between the 
Miras in the two populations. Especially in the log P-(J — H)o dia- 
gram, the difference between the instruments seems to have a larger 
effect than any possible population differences. We will therefore 
use the PLR, eqs.QJ{3] obtained from the observations by 104. The 
possible population effect will be discussed again in section |531 



5.2 Estimation of fj,o and Ak s 

We can get the apparent distance modulus for an individual Mira 
of known period using the PLR (eqs.[T}[3j if a mean magnitude is 
available in a given filter. The apparent distance moduli in JHK S 
are related to the true distance modulus /10 and the extinction Ak s 
as follows: 



Ho + rjA Ka =J-Mj, 
(io+r H A Ks = H-M H , 
fio + A Ks =K B - M K , , 



(4) 
(5) 
(6) 
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Figure 19. Period-colour relations for Miras obtained in earlier 
papers. Line st yles and wid t hs for differ ent samples are ind icated 
at the top: F894Feast et alj |l989t). G95-lGlass et alj Jl995l) . 104- 
llta et alj J2004bf ). and M05- lMatsunaga et alj J2005h . The projec- 
tion of each line on the horizontal axis indicates the period range of 
the samples used to derive the relation. 



where r.j and th are defined as Aj/Ak b and Ah/Ak b re- 
spectively. The extinct ion coefficients have been obtained by 
Nishivama ~et~aLl d2006t^ >. 



rj = 3.021, r H = 1.731, 



(7) 



using the same instrument as ours. When a period and mean magni- 
tudes in two or more filters are available, we can solve for the two 
unknown quantities (/io, Ak s )- For a Mira with all three JHK B 
magnitudes available, we can make three estimates of (/io, Ak b ) 
by using the three different pairs of filters. In order to specify the 
band pair used to derive a particular (/io, Ak b ), we use the band 
names as superscripts, e.g. /j,o IK and Ak s k . For example, the so- 
lution of eqs. (O and ([6} can be written as follows: 

H — K s Mb- M Kb 



A, 



rn - 1 rH — 1 
r H K s - H r H M Ks - M H 



1 



(8) 
(9) 



We have made distance estimations using only Miras with 
P ^ 350 d, amplitudes not smaller than 0.4 mag and at least 
two mean magnitudes having Mflag = 0. There are 175 Miras 
fulfilling these requirements; the numbers for the three pairs of fil- 
ters, being 143 for (H, K s ), 37 for (J, K s ) and 69 for (J, H). Ta- 
ble [9] lists the (/io, Ak s ) values found. For the 37 Miras which 
have all three JHK S means, we were able to obtain three esti- 
mates of (/io, Ak b ). Ideally, the values found from the different 
pairs should agree with each other. The (/io, Ak s ) values are mu- 
tually compared in Fig. [20] The points are distributed around the 
diagonals but show some scatter. 

In the following, we estimate the uncertainties in solving 
(/io, Ak s ) and show that the scattering in Fig.|20]can be explained 
quantitatively. The solutions ©, and similar ones for other 
pairs of bands are affected by the uncertainties in our photometry 



Table 9. Extinctions Ak and distance moduli /io estimated for Miras. 
Three estimates from different pairs of filter bands are listed whenever 
available. We put 99.99 if they are not available. The full version will 
be available in the online journal. 



No. 


A HK 




A JK 




A JH 




0002 


3.19 


14.22 


99.99 


99.99 


99.99 


99.99 


0004 


99.99 


99.99 


99.99 


99.99 


1.52 


14.40 


0010 


2.69 


14.35 


2.70 


14.34 


2.71 


14.32 


0016 


4.47 


14.07 


99.99 


99.99 


99.99 


99.99 


0023 


1.71 


14.96 


1.82 


14.84 


1.89 


14.65 


0034 


4.05 


14.30 


99.99 


99.99 


99.99 


99.99 


0035 


2.69 


14.87 


99.99 


99.99 


99.99 


99.99 


0039 


99.99 


99.99 


99.99 


99.99 


1.59 


14.68 


0040 


2.02 


14.61 


2.03 


14.61 


2.03 


14.60 


0074 


99.99 


99.99 


99.99 


99.99 


2.00 


14.46 



(JHK S ) and in the predictions from the PLR (Mj,M H ,M Ks )- 
The prediction of absolute magnitudes based on the PLR cannot 
be totally accurate because there is scatter around each relation as 
seen in Fig.[T8] The scatter can be as large as 0. 19 mag, but the de- 
viations are not mutually independent. For the 134 LMC Miras in 
Fig. [18] the deviations from the PLR in J, H and K s are strongly 
correlated, which should be borne in mind when estimating the 
uncertainty caused by the width of the PLR. The uncertainties in 
the JHK S magnitudes are considered to be mutually independent, 
and we estimate that the errors in the mean magnitudes of JHK S 
are 0.02 mag. We calculated expected uncertainties in the differ- 
ences between estimates of Ak s or /io, such as A^ — A^ and 
/i ,JH — Hq JK . The results are listed in Table 1 101 together with the 
scatters actually obtained in our results. The predicted uncertainties 
agree with the actual values well. The uncertainties in the estimated 
periods can also affect the estimation of (/io, Ak s )- However, they 
produce errors of almost the same size for different band pairs, and 
they compensate each other when we take the differences listed in 
Table 1 101 The extinction coefficients rj and th averaged over the 
GC region are reported to have uncertainties of about 1 percent 
dNishivama et al.|[2006al) . From eqs. ([Sj and l[9), the resultant er- 
rors in A™ and/i Hif are 0.06 and 0.07 mag, respectively. How- 
ever, they do not affect the differences between the results from 
each colour combination. We therefore conclude that the scatters 
in Fig. [20] are adequately explained by the expected errors in our 

met hod. 

iGlass et al. suggested that the Miras in the Sgr I field 

have larger (H — K)o and smaller (J — H)o than those in the 
LMC, at a given period, as can be seen in Fig. [19] If such a dif- 
ference exists, our estimates would give larger A^ K and s maller 
A^ values for the bulge Miras. In effect, in IGlass et all d 1995b 
the estimated A^ K values are systematically larger than the Aj^ 1 



by 0. 1-0.2 mag, depending on the period range considered. We did 
not, however, find such a shift in Fig. [20] It is therefore suggested 
that the PCRs of the Miras in the LMC and the bulge are similar to 
each other in the IRSF/SIRIUS system. This may be confirmed in 
the future from IRSF photometry of the Miras in the Sgr I field and 
NGC 6522. 

We have plotted histograms of Ak s in Fig. [2T] and those of 
/io in Fig. [22] where each panel corresponds to the result for each 
pair of filter bands. An obvious difference between the results from 
different pairs is that Miras with Ak s > 2.5 mag are well detected 
in panel (a) of Fig.[2T]but not in panels (b) and (c). This is because 





14 15 

HK 

Figure 20. In the above three panels, we compare Ak values obtained with magnitudes in different pairs of filters. Comparisons for values 
are also plotted in the lower panels. 



Table 10. Errors in some values related to the estimation of 
(fj.0, Ak b )- The predicted uncertainties are based on eqs. (8)-{9) 
and similar solutions. The parts related to observational magnitudes 
( JHK B ) and those related to the PLR predictions are individually 
listed together with the totals. The scatters found in Fig. 1201 are 
listed in the last column. The uncertainties in the periods affect the 
values given in the last two lines, but they are not included here. 
They have almost no effect on the other six lines (see text). 



Value 



Predicted uncertainty Observed 
Obs. PLR Total scatter 



A JH 

A K. 

A " K 

A K. 



A JK 

K f 
A -' H 



A JK — A 



0.019 
0.053 
0.034 
0.038 
0.092 
0.034 



0.038 
0.105 
0.067 
0.111 
0.182 
0.067 



0.042 
0.118 
0.075 
0.117 
0.203 
0.075 



0.040 
0.110 
0.070 
0.123 
0.192 
0.070 



A l K 



0.039 0.082 0.091 
0.039 0.179 0.183 



such highly-reddened Miras are too faint in J for our survey. In fact, 
none of the Miras with A^^ K > 3 mag have J-band magnitudes. 




5.2.7 A note on small-amplitude variables 

We did not take periodic variables with small-amplitudes (A J < 
0.4, AH < 0.4, or Aif s < 0.4) to be Miras. Among the sources 
with periods between 100 and 350 d, 27 are of small- amplitude 
in H and/or K s , i.e. AH < 0.4 and/or AK S < 0.4. If we esti- 
mate their distances assuming that all of them are Miras, the re- 
sultant distribution of fiQ 11 ^ looks different from that in the panel 
(a) of Fig. [22] In the case of Fig. [22] there are almost no objects at 



Figure 21. Histograms of Ak b values obtained by using magni- 
tudes in different pairs of filters. 



Hq IK < 14. In contrast, five out of the 27 small- amplitude vari- 
ables would be considered to be closer than (J-o IK = 14 mag using 
(|2j, ([5J and J9). Their distance moduli could be around 14.5 mag 
or slightly larger than that if they foll ow the sequence B or C' in 
the log P-K s diagram jlta et al.l20 04a) rather than the sequence C. 
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Figure 22. Same as Fig.HT] but for no- 
There is no reason to suppose that small-amplitude Miras are pref- 
erentially located closer to us; thus it is likely that at least these five 
objects are not Miras. We decided for safety to exclude the small- 
amplitude variables from the list of Miras (even though many of 
them may follow the sequence C). As it happens, our estimation of 
the distance to the GC is not affected even if we include the small- 
amplitude variables in the following discussion. 



5.3 Errors in estimation of (fio , Ak s ) 

As has been discussed in the preceding section, the solutions of {8} 
and ([9} are affected by the uncertainties in our photometric error 



and the width of the PLR. We calculated their effects on A I 



and 



H and listed them in Table 1 101 We now summarize the sources 
of error in estimating (fj,o, Ak s ), when both of the above are in- 
cluded, in Table 1111 We classify the error sources into three types: 
random, systematic and statistical. The first type relates to individ- 
ual Miras, while the second has a common effect on all the Miras 
in our sample. When we study the distance to the GC, the statisti- 
cal error carries more weight than the random errors for individual 
objects because we make use of many objects. 

Our estimated periods have uncertainties of 0.025 dex in 
log P as seen in section [3T2l They introduce uncertainties of 0.07- 
0.09 mag in the predicted absolute magnitudes via eqs. {TJ— <[3j. 
Since errors added to Mj, Mh, and Mk b have the same sign, 
their effect on A^ b k is almost negligible (~ 0.01 mag). On the 
other hand, Ho HK is affected by ±0.09 mag. This error is random 
and must be included with the photometric error and that from the 
width of the PLR. The result is a total error of 0.21 mag which 
randomly affects the distance estimates for individual Miras. 

There are also uncertainties i n the zero points of both the 104 
catalogue and the photometry of iNishivama et all J200 6a) which 
we used for the calibration. The difference between the zero points 
in these works is expected to be small because both were cali- 
brated against standard stars in lPersson et al.l dl998t> . The standard 



star #9172 was used by INishivama et all d2006al) and the uncer- 
tainty of their m agnitude scale is rather small (see section 2.1 in 
INishivama et alj |2009). On the other hand, 104 used a few dozen 
of Persson's standards. They used A'-band, not I"f s -band, magni- 
tudes of these standard stars, but the differences are negligible. The 
zero-point uncertai nties were not given, but are not expected to be 
larger than those in lNishiyama et all j2006al) . We therefore adopted 
0.01 mag for the uncertainties in the magnitude scale. 

The extinction coefficients r j and th give rise to uncertainties 
of 0.06 mag in 

Ak K and 0.07 mag in Hq HK as already mentioned. 
Here we do not consider any variation of extincti on coefficient 
with th e line of sight in our survey region. Recently, Gosli ntTet al.l 
d2009l) have suggested that there are substantial variation on a very 
small scale (~ 5 arcsec). Such a variation would lead to different 
values of Ak s from different pairs of filters. In contrast, the scatters 
in Fig.|20]can be explained without the variation. This suggests that 
r j and rg are almost constant in our region. 

We need at this point to consider the uncertainties arising from 
the calibration of the PLR. We derived eqs. (T|l-l[3) assuming the 
distance modulus to the LMC to be ^lmc = 18.45 mag. The 
A*lmc value is one of the most fundamental values in the cosmic 
distance scale and has been widely discussed. Most of the recent 
measurements are around 18 .4-18.5 mag (e.g. Ivan Leeuwen et al.l 
120071 ; iGrocholski e7aT1l2007l : Clement, Xu & Muzzin, 2008). We 
adopt ^lmc and its systematic error as being 18.45±0.05 mag. A 
recent result based on the PLR of type II Cepheids in the LMC 
is also consistent with this value (Matsunaga, Feast & Menzies, 
2009). Note that this error has nothing to do with the estimation of 
Ak b . 

It may be asked whether there is a population effect on 
the PLR of Miras tha t might also introduce a systematic error. 
IWhitelock et al.l (2008) has discussed the present status of the cali- 
bration of the PLR although they focus on the log P-K relation and 
do not include other wavelengths. They compared the PLR of the 
LMC Miras with those in the solar neighbourhood, globular clus- 
ters and the Galactic bulge. The relations for these environments 
agree with each other within an uncertainty of around 0.1 mag. 
Unfortunately, this uncertainty includes the error in /^lmc because 
one cannot make the above comparisons without its use. Their con- 
clusion is that population effects on the PLR are negligible within 
the current accuracy. As for the period-colour relation, some au- 
thors have reported systematic differe nces between Miras in dif- 
ferent stellar populations. For example. iFeastl ( fl996l) compared the 
log P-( J — K)o relation of Miras in different stellar systems and 
found a difference as large as 0.05 mag in ( J — K)o in the period 
range between 100 and 350 d. A shift of 0.05 in ( J — K)o would 
lead to a shift of 0.025 mag in the estimate of (/io, Ak s )- We will 
assume here that the population difference between the LMC and 
the Galactic bulge may introduce a colour difference of 0.05 mag 
mH — K s based on the top panel of Fig.[l9]and that the zero point 
of the PLR in K s is not affected. Thus uncertainties of 0.07 mag 
are introduced into both A^ s K and Ho HK . 

In summary, there are four major sources of systematic error 
as listed in Tablel 1 1 1 and their total size is expected to be 0. 1 1 mag. 

5.4 Distance to the GC 

We have discussed the uncertainties in our distance estimation in 
the preceding section. While random errors in /io were found to be 
about 0.2 mag, the peak at around no = 14.5 mag has a standard 
deviation of 0.19 mag, which is obtained for the objects between 
14 and 15 mag in Fig. [22] Therefore, the width of the distribution 
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Table 11. Sources of error and their amounts for AS K and 



Error source 



Error size (mag) 



A HK 



Random uncertainty 



photometric error 


0.04 


0.04 


width of the PLR 


0.08 


0.18 


period 


0.01 


0.09 


total 


0.09 


0.21 


Systematic uncertainty 




photometric zero point 


0.01 


0.01 


extinction coefficient 


0.06 


0.07 


the LMC distance 


0.00 


0.05 


population effect on the PLR 


0.07 


0.07 


total 


0.09 


0.11 



Statistical uncertainty 
standard deviation of mean 
total 



0.02 
0.02 



can be explained by the random errors. This carries the implication 
that the Miras are strongly concentrated to the GC. 

Our survey field extends to approximately 12 arcmin in ra- 
dius, corresponding to 30 pc at the distance of the GC. This region 
is well within the so-called Nuclear Bulge which is distingui shable 
from the large-scale Galactic bulge (Serabyn & Morris 199rj). The 
Nuclear Bulge consists of an R~ 2 stellar cluster which contains 
a significant stellar population as young as 10 7-8 yr. Long-period 
Miras (say P > 600 d) and OH/IR stars indeed show a strong con- 
centration towards the central 100 pc region (Lindqvist, Habing & 
Winnberg, 1992; Wood, Habing & McGregor, 1998). However, it 
is known that short-period Miras are relatively old (see the review 
bv lFeastll2009h . so that the Miras with P = 100-350 d whose dis- 
tances we obtained may have a different spatial distribution from 
the younger populations. Frogel, Tiede & Kuchinski (1999) showed 
that the younger component of the stellar population is more con- 
centrated than the overall bulge population. Nevertheless, a large 
fraction of the old stellar population is also located in a rather nar- 
row region. As a rough estimate, more than eighty percent of ob- 
jects projected in our fields of view are expected to be within 300 pc 
of the GC if we assume that their density distribution follows the 
truncated power law used in Binney, Gerhard & Spergel (1997, see 
their eq. 1). In contrast, the standard deviation around the peak in 
Fig. 1221 corresponds to ±700 pc. We therefore conclude that most 
of the Miras that we have detected lie at the same distance and are 
centred around the GC within our accuracy, so that the distance to 
the GC can be derived legitimately from their distance statistics. 

It is necessary to consider any further biases which may have 
arisen. As mentioned in section [3~3l we could not always detect the 
Miras in all three filters. The distances which we can estimate are 
limited by the observable magnitude range and are also by the ex- 
tinction Ak s ■ Based on the limits in JHK S (Fig. [3), it is possible to 
define the range of (fio, Ak b ) over which one can observe objects 
with a certain absolute magnitude. We here adopt Mj — —5.0, 
Mh = —5.6 and Mk = —5.9 as typical magnitudes for Miras 
with P — 100-350 d. In Fig. [23] we plot boundaries which cor- 
respond to the detection limits (dashed curves) and the saturation 
limits (solid curves) in the (/in, Ak s ) plane. The boundaries are 
different for different filters as shown in Fig. [23] The right-hand 
panel presents the same curves as the left panel but with a linear 
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Figure 24. Distribution of the n HK values obtained. The error 
bar indicates the size of the Poisson noise (v~N) for each bin. The 
mean value and width of the fitted Gaussian are indicated in the 
panel. 



distance scale to show more clearly the range we are interested in. 
The J-band detection limit implies that we cannot detect Miras if 
they are more heavily extincted than 2-2.5 mag in Ak b ■ Note that it 
is possible to detect Miras over a wider range of distance even in the 
J-band if the reddening is smaller than Ak b = 2 mag. This is why 
the differences between the distributions are not large in Fig. [22] 
One can observe the Miras with larger Ak s by using H and K B , 
and so the number of (/j, , A K ) estimates is significantly larger 
than when J is involved. We therefore use only fi Q HK andAg K to 
discuss the distribution of Miras in what follows. The shaded areas 
in Fig. 1231 indicates the range of (fio, Ak s ) where we can detect 
Miras in both H and K s . When we investigate the distance distri- 
bution of Miras, we exclude those with Ak b > 3 mag because the 
Miras on the far side of the Galactic bulge may not be detected at 
such high extinctions. One should bear in mind that such a limita- 
tion can introduce a selection bias in the values of (/Lto, Ak s ) that 
we found. The rectangle in the right-hand panel indicates the range 
we made use of in determining the distance to the GC. 

We may now proceed to estimate the distance to the GC 
based on the 100 Miras which have estimated values of [1q IK be- 
tween 14 and 15 mag and Ar between 1 and 3 mag. The lat- 
ter condition is set to avoid selection bias in the (fio, Ak b ) plane 
(Fig.|23t. One can expect the spread in the estimated distance mod- 
uli to accord with the random errors. As shown in Fig. [24] the 
distribution is well fitted by a Gaussian function with a peak of 
fio = 14.58 mag. The peak and the width from the fitting are in- 
sensitive to the bin width. As already mentioned, random errors 
for individual Miras are not important when we discuss the dis- 
tance to the GC because they average out. We take the standard 
deviation of the mean, viz. 0.02 mag, as the uncertainty in the dis- 
tance. As a result, the distance modulus to the GC is estimated to be 
14.58±0.02 (stat.) ±0.11 (syst.) mag where the systematic error 
is included ( section [53V We denote the linear distance as Ro and 
compare our result, Rq = 8.24±0.08±0.42 kpc, with those of oth- 
ers. When necessary, we re-calculate Ro by using /ilmc = 18.45 
(if the distance scale in a given work depends on this quantity). 

In earlier works Ro values have been obtained by various 
methods. That involving stellar orbits around the central black hole 
Sgr A* is a potentially powerfu l tool, but there still re mains a 
significant uncertainty. Al though lEisenhauer et alj d2005h derived 
Ro = 7.62 ± 0.35 kpc. iTrippe et alj d2008h has found that sys- 
tem atic errors neglecte d in the 2005 paper may change the re- 
sult. iGhez et alj j2008h obtained Ro = 8.0 ± 0.6 kpc using the 
same method but based on a different set of data. The same au- 
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Figure 23. Schematic diagrams to show the parameter space of (fio, Ak s ) where we can detect Miras in JHK B . Dashed 
curves show the detection limits while solid curves show the saturation limits for JHK B as indicated. In the online version, 
blue, green and red curves corresponds to the limits in J, H and K B respectively. The two panels plot the same curves, but 
the y axis in the right one is on a linear scale. The rectangle in the right-hand panel indicates the region we use to estimate 
the distance to the GC. 



thors obtained R = 8.07 ± 0.32 (stat.) ± 0.13 (syst.) kpc as a 
statistical parallax using the kinematic properties of several hun- 
dred stars around the G C. Results from other methods include: 
Nishiyamaet al. (2006b; Ro — 7.52 kpc, red clump stars), Groe- 
newegen, Udalski & Bono (2008; Ro — 7.75 kpc, type II Cepheids 
and RR Lyr variabl es). Our result is simil ar to the recent results by 
iTrippe et all J2008h and lGhez et al] J2008I) . 

A few other auth ors have estimated R p from Miras. Pioneer- 
ing work was done by [Glass & Fe ast ( 1982), who conducted near- 
IR photometry for a few doze n of Miras t oward s the Baade win- 
dows (Sgr I and NGC 6522). iGlass et al.l dl995h investigated the 
Miras towards the Sgr I field after further monitoring. Their result 
was Rq — 8.3 ± 0.7 kpc, but they mentioned that their Ro value 
will increase by 0.2 kpc if the Sgr I field is part of b ar at ~ 45° to 
the line of sight. Groenewegen & Blommaert] d2005h used the Mi- 
ras found in the OGLE-II survey fields which cover various lines of 
sight between ±11° in I, mainly at b ~ —3.5°. They found that the 
difference between the distance moduli of the LMC and the Galac- 
tic bulge is 3.72 mag, which corresponds to Ro = 8.9 ± 0.7 kpc. 
They also argued that Ro can be as small as 8.35 kp c if one as- 
sumes a strong metal dependency of the PLR of Miras. IGlass et al.l 
(2009) adopted R = 8.5 kpc by comparing their plots of mid-IR 
magnitudes against periods for pulsating red giants in the LMC and 
the Galactic bulge including semi-regulars on other sequences. The 
previous results show rough agreements with ours considering the 
systematic uncertainties. Note that the fields they examined are at 
least a few degrees from the GC. Our report is the first to estimate 
Ro by using the Miras towards the GC region itself. Furthermore, 
the method described here has made it possible to overcome the 
complicated effects of the interstellar extinction. 



5.5 Interstellar extinction 

In the course of this work, besides their distances, interstellar ex- 
tinctions have been estimated for more than 150 Miras. Because we 
estimate both at the same time, this method has the potential to in- 
vestigate the extinction in three dimensions directly, revealing how 
it changes with distance and with line of sight. In spite of the fact 
that our objects are located at the same distance, the Ak s values 
cover a large range; the foreground extinction varies strongly with 
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Figure 26. A £ K are plotted against fi H 



the line of sight. In Fig. [25] we plot the locations of the Miras with 
values in galactic coordinates. The sizes of the symbols vary 
according to Ak 3 as indicated. Whilst circles with various sizes 
scatter across our observed area, there are some clumps of similar- 
sized circles. 

Fig. I26I plots against fJ.p HK ■ The cross indicates the 

sizes of the errors in the two quantities. Although most of the sam- 
ples in this study are located at the same distance within the accu- 
racy, there are two sources closer and two further than the main 
group. The Ak b values obtained for these objects do not show 
systematic differences from those for the objects with no = 14- 
15 mag. Unfortunately, the number of tracers with (fio , Ak s ) is too 
small to carry out a detailed analysis of the three-dimensional ex- 
tinction map. Also note that our sample is limited in the (fio , Ak s ) 
plane (see Fig. I23t. Those with smaller extinction and closer to the 
Sun will be saturated in our survey, while those with larger extinc- 
tion and further than the GC are expected to be too faint. 

We here estimate the foreground interstellar extinction to- 
wards Sgr A* (I = -0.05575°, b = -0.04617°). Fig. [27] shows 



the locations and the A^ s K values (and fi H K in parentheses) for 
the Miras around Sgr A* . Four Miras are located within a radius 
of 1.5 arcmin, indicated by the dashed circle. It is suggested that 
the Ak s in the direction towards Sgr A* is approximately 3 mag. 
Our value corresponds to Ay ~ 45 mag if we use the extinction 
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Figure 25. Locations of 143 Miras with A^ K obtained. The obtained extinction for each object is indicated by the size of a circle 
as shown besides the panel. 



coefficient of Av /Ak„ = 16 obtained bv lNishivama et al] J2008h . 
In contrast, extinction v alues of 2.6 and 2.8 ma g in Ak„ were ob- 
tained toward Sgr A* by Eisenha uer et al.1 feOQSl) and lSchodel et al.1 
J2007h with higher angular resolution. The latter found that the ex- 
tinction changes on rather small scales. Our estimate should there- 
fore be regarded as a typical value for a larger region surrounding 
Sgr A* rather than the value in its precise direction. 



5.6 Mid-IR PLR 



iGlass et alj d2009h recently investigated the PLR of Miras in the 
mid-IR range. Their samples consist o f the Miras towards the LMC 
and the NGC 6522 Baade window. ISchultheis etall < |2009|) fitted 
the mid-IR PLRs for the G01 Miras around the GC, but there re- 



main large scatters because it was not possible to correct for the 

variable interstellar extinction. 

W e have plotted the IRAC magnitudes from iRamfrez et al.1 
(2008) against the periods in Fig.|28|for our Miras with (jUn, Ak^)- 
The fitted linear relations are: 

[3.6] = -4.14(±0.20)(logP- 2.3) - 7.32(±0.03), (10) 
[4.5] = -4.03(±0.20)(logP - 2.3) - 7.48(±0.03), (11) 
[5.8] = -4.33(±0.24)(log P - 2.3) - 7.88(±0.04), (12) 
[8.0] = -4.46(±0.33)(log P - 2.3) - 8.14(±0.06), (13) 

indicated by the solid lines. The residual standard deviations around 
( fTOt . ( flTi and ([T3ll are 0.29, 0.29, 0.32 and 0.42 mag, respec- 
tively. We here used the (1q IK and A^ b k obtain ed for individual 
Miras and the extinction coefficients obtained by [Nishivama et all 
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Figure 27. This is similar to Fig. [25] but is a closeup for the prox- 
imity of Sgr A* which is indicated by the cross. The dashed circle 
has a radius of 1.5 arcmin and a centre at Sgr A*. The Ajf s (and 
/jo in parentheses) is indicated for each object. 
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Figure 28. Period-luminosity relation in SpitzerflRAC filters. 
The fitted r elation s are shown as solid lines, while the PLR in 
iGlass et al] j2009h are indicated by dashed lines. Note that they 
adopted a distance differing from ours by 0.12 mag and we have 
shifted their fits by this amount. 



< l2009h to get their absolute magnitudes. Because most of the Miras 
are found at the same distance within our accuracy, it is also rea- 
sonable to use the common distance modulus of 14.58 mag. Even if 
we do so, the resultant relations do not change from eqs. dlOM13t 
including their scatters. 

The scatter around the PLRs that we have found are rather 
small considering the uncertainty in the (no, Ak s ) and the fact 
that the IRAC magnitudes are based on a single-epoch survey. Un- 
certainties in the apparent distance moduli in the IRAC filters are 
estimated at 0.18 mag for individual Miras. The mid-IR amplitudes 
of Miras are not well kno wn, but those at a round 3.5 fim can be 
as large as in the .ft. -band dFeast etal.ll 19821) . These uncertainties 
suggest that the intrinsic widths of the mid-IR PLRs are narrower 
than in Fig. [28] 

In Fig. [28] we also plot the relation obtained for the C- 



sequence bv lGlass et al] [2009) as the dashed lines. Their adopted 
distance, 14.7 mag for the NGC 6522 field, differs from ours by 
0.12 mag and was compensated for before plotting. The slopes 
agree well with ours. In contrast, the zero points are slightly shifted 
depending on the wavelength; approximately by +0.3 mag for the 
log P-[3.6] relation and —0.2 mag for the log P-[8.0] r elation. The 
reason for the different zero points is unclear. Note that lGlass et al] 
J2009I) included semi-regulars and long-period Miras which may 
have mid-IR excesses in their sample. 



5.7 Future prospects 

As has been shown, it is possible to obtain both distances and ex- 
tinctions for Miras with P — 100-350 d by making use of the 
PLR in two or more near-IR filter bands. This is a very promising 
tool for investigating Galactic structure. One of the major sources 
of error is the residual scatter around the PLR. Because we do not 
have a priori information on how a certain Mira deviates from the 
fitted PLR, the scatter introduces an uncertainty in the distance es- 
timation. We expect that the intrinsic scatter is smaller than what 
we found in Fig. [T8] The magnitudes used in Fig. [T8] are intensity 
means of ten-epoch monitoring data, which can introduce an un- 
certainty as large as 0.1 mag based on the Monte-Carlo simulation 
presented in section [33] The scatter should get smaller if we take 
mean magnitudes more accurately. 

The uncertainty in the zero point of the PLR and the pos- 
sible population effect may also introduce significant errors; it is 
important to note however that these will be systematic ones. Fu- 
ture parallax observations of Miras in the solar neighbourhood ma; 



hood may 
ail d2003h 



red uce these uncertainties. For example. [ Vlemmings et 
and IVlemmings & van LanHeveld3 d2007t) have reported the par- 
allaxes for five Miras with OH maser emission. As reviewed by 
IWhitelock et al] d2008l) . their results are consistent with the PLR 
of the LMC Miras within their uncertainties, which are still not 
small enough to draw strong conclusions about the zero point or 
the population effect. VERA, a Japanese VLBI project, is one of 
the most promising projects for calibrat ing the PLR of Miras in 
the Galaxy (Nakagawa et al. 20081 120091) . However, the use of ra- 
dio parallaxes is constrained by the small numbers of Miras with 
maser emission in the short-period range, say P ^ 300 d. In our 
method for deriving the distance and the extinction, we make use 
of short-period Miras to avoid the effects of circumstellar extinc- 
tion and/or hot-bottom-burning. In this reg ard, future optical and 
near-IR ast rometric satellites s uch as GAIA dTuron et al.l2008h and 
JASMINE dGouda et al]|2008l) can play a fundamental role in the 
calibration of the PLR. 



6 SUMMARY 

This paper reports 1364 variable stars in a 20' x 30' region around 
the GC found through our near-IR survey covering seven years. We 
have obtained periods for 549 objects, most of which are considered 
to be Miras. We have described a method for obtaining distances 
and extinctions for Miras with periods between 100 and 350 d by 
making use of their PLR. The technique depends on knowing mean 
magnitudes in two or more filters in the near IR. A particular effort 
has been made to assess the uncertainties; the random error is about 
0.2 mag in the distance modulus for an individual Mira, while the 
systematic error is about 0. 1 1 mag, mainly originating from a com- 
bination of uncertainties in the LMC distance (±0.05 mag), the 
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extinction coefficient (±0.07 mag) and the population effect on the 
PLR (±0.07 mag). 

We obtained the distances of 143 Miras by using the H- and 
its-band mean magnitudes. They are strongly concentrated to the 
GC, whose distance modulus was estimated at 14.58 mag. The sta- 
tistical error is 0.02 mag from the standard deviation of the mean 
and the systematic error is again 0. 1 1 mag. Our statistical error is 
the smallest among estimates to date, and the distance we have 
obtained shows rough agreement with recent results from other 
methods, within the systematic uncertainties. The PLR of Miras 
is shown to be a powerful tool for studying Galactic structure, al- 
though there remains a significant systematic error. It is expected 
that better calibration of the PLR will improve the accuracy of the 
method; this may be achieved by parallax projects such as VERA 
and GAIA. 

We also have discussed the distribution of the interstellar ex- 
tinction towards the GC. Because we can estimate distances and ex- 
tinctions for individual Miras, they can be used as tracers to reveal 
three-dimensional structure. It is possible to study how the extinc- 
tion changes with distance along various lines of sight, although 
the low space density of Miras may prevent us from investigating 
small structures. We have confirmed that the extinction towards the 
GC region varies complicatedly between 1.5 and 4 mag in Ak b , 
excepting towards extremely dark nebulae, and we have estimated 
the extinction towards the Sgr A* region to be Ak s ~ 3 mag. 
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